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This is an extended Reply to Comment by A. Sergeev, M.Y. Reizer, and V. Mitin jarXiv:0906.2389| 
on our Letter [Phys. Rev. Lett. 102, 067001 (2009)]. We explicitly demonstrate that all claims 
by Sergeev et al. are completely unfounded, because their underlying theoretical work contains 
multiple errors and inconsistencies. For this reason, there is no need to revise the existing theories 
of thermoelectric response in superconductors. 



PACS numbers: 74.40.+k, 72.15.Jf, 74.25. Fy 

INTRODUCTION 

In a series of recent papers (see, Ref. Q and references 
therein), Sergeev, Reizer, and Mitin (SRM) have ar- 
gued that all existing results on thermoelectric response 
in fluctuating superconductors are fundamentally flawed 
and must be revisited. According to SRM, these pre- 
vious allegedly incorrect works include the pioneering 
works of UUah and Dorsey ^] and those of Ussishkin, 
Sondhi, and Huse 0-01) the relevant chapter in the book 
of Varlamov and Larkin our recent Letter_0|, and 
the related work of Michaeli and Finkelstein [Sj. Note 
that these researches have employed at least four 

different approaches to calculate the Nernst effect in a 
superconductor: Time-dependent Ginzburg-Landau the- 
ory [l,!!,!!!, microscopic techniques based on the Matsub- 
ara diagrammatic technique [J 0] and on the quantum 
kinetic equation and finally simple phenomenologi- 
cal arguments, which relate the Nernst coefficient to the 
temperature dependence of the chemical potential for the 
carriers ■ All these approaches consistently provide the 
same order of magnitude for the Nernst coefficient, which 
has been shown to much exceed that in a Fermi liquid 
(hence, we called the effect "giant" in our Letter). Such a 
giant Nernst signal has been detected in a variety of well- 
known experiments in the high-temperature cuprates 
and conventional superconducting films l3, llj- There- 
fore the existence of giant Nernst effect constitutes a solid 
experimental fact, which has been the actual motivation 
for the aforementioned theoretical works. However, SRM 
"strongly object" to the existence of these results and ex- 
perimental facts, arguing that the "numerous recent 
theories grossly overestimate the thermomagnetic coeffi- 
cients" and calling them in the abstract of their recent 
comment [l[ "ridiculously large." These strong claims of 
Sergeev et al. are based solely on their own alternative 
calculation approach, sketched by them in Ref. 12 1. 

In 2009 SRM had escalated the concern about the ex- 
isting theories of the Nernst effect, by posting a comment 



[H on our Letter. Despite SRM's evident errors, their 
Comment has recently been accepted for publication in 
Physical Review Letters. This manuscript is an extended 
version of our Reply. 

Our Reply contains the following: (i) In the first part, 
we analyze the paper [12], which is needed because all 
SRM's criticism on the existing theories is based on this 
single paper. Hence, a careful analysis of Ref. 12| is the 
only means for us to refute SRM's criticism of the works 
by us and others. As a result of this exercise, we are 



able to show in the first part of the Reply that Ref. |12| 
contain serious errors. We identify the main problems 
of SRM's treatment as likely originating from the combi- 
nation of inconsistent use of a specific gauge in a single 
calculation, missing pieces in the relevant diagrams, and 
technical mistakes, (ii) Next, we suggest specific tech- 
nical steps to remedy these problems of SRM and bring 
their approach in accordance with the exiting theories. 
But most importantly, we conclude that since the paper 
12l | is manifestly incorrect, there is no need to revise all 
previous existing theories, (iii) In the last part, we elab- 
orate on the phenomenological Eq. (1) that we suggested 
in our Letter Q as a simple intuitive argument behind 
the experimentally observed giant Nernst effect. 



ANALYSIS OF PHYS. REV. B 77, 064501 (2008) 
BY SERGEEV ET AL. 



Ref. 



12| concentrates entirely on the discussion of "a 
gauge-invariant microscopic approach," but we show be- 
low that it is the very gauge-invariance that is explicitly 
violated in the SRM's treatment. Below we choose to 
focus on two specific technical problems (among many) 
that suffice to prove the inadequacy of their proposed 
calculational method and therefore make the criticism 
mute. 
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On the calculations of the Nernst effect in a normal 
metal by Sergeev et al. 

One of only a few existing and commonly-accegted re- 
sults on the Nernst effect is that of Sondheimer [ISj , who 
first calculated the effect in a Fermi gas back in 1948. 
Any technique that claims to be of relevance to more 
complicated Nernst phenomena, such as those due to su- 
perconducting fluctuations, must recover Sondheimer's 
formula as a basic "sanity check". In Sec. II of their 
work, SRM do consider the case of a normal metal but 
fail to explain the origin of Sondheimer's [ist formula in 
a comprehensible way. In addition, their discussion con- 
tains an error, which propagates into further analysis. 

For noninteracting fermions, the thermoelectric re- 
sponse is described by the Kubo-like diagram constructed 
of the heat and electric vertices and two Green functions. 
In the absence of a magnetic field, the heat vertex is 
^pV = ^pp/m. For nonzero magnetic field, both the ki- 
netic energy ^p and the momentum p itself, should be 
modified by the vector potential A. SRM forget about 
the latter and thus come to a wrong expression for the 
heat current (Eq. (5) of Ref [l2|): 



'SRM — 

P 



v^pa+ap + ^— (v- A)a+ap. (1) 



Further they forget to incorporate the vector potential 
A into the electric current to maintain gauge invariance. 
Then the correlator {j'^j'^) is claimed to be calculated, 
but (i) a gauge in which calculations are performed is 
poorly defined and (ii) the external momentum is unspec- 
ified (below, we explain why it is important). We reiter- 
ate here that while SRM concentrate their work as well 
as their criticism of others on the gauge-invariance issue, 
which indeed is important, they fail to satisfy the gauge 
invariance in their own work, even at the non-interacting 
level. Their "gauge-invariant" expressions contain ex- 
plicitly a vector potential, which is inserted into selected 
pieces of diagrams in an uncontrolled fashion and calcu- 
lations are performed in an unspecified gauge. 



On the calculations of the Nernst effect in a 
fluctuating superconductor by Sergeev et al. 

In Sec. II of Ref [l2| Sergeev et al. analyze the heat 
current transferred by fluctuating Cooper pairs in the 
vicinity of Tc. They correctly obtain that the heat cur- 
rent is proportional to the gauge-invariant momentum 
(Eq. (11) of Ref. [i3): 



B^RM oc w(q + 2eA/c) 



(2) 



and they claim that it is the term with A that is their key 
new finding and that this term was allegedly overlooked 
in all previous calculations. However, SRM again fail to 



include A in the electric vertex and draw the diagrams, 
extracting A from the propagators and from the heat ver- 
tex ([2]) only. Then they say: "Further calculations of the 
diagrams for Cooper pairs are similar to that for nonin- 
teracting electrons", thus inheriting all the inconsisten- 
cies from their normal-metal treatment. Even without 
that, the expressions determining current ([T]) and ([2|) are 
different in their vector structure: the first one contains 
v(v • A), while the second does not! 



Correcting the problem with gauge invariance 

Here we elaborate on the main technical problems of 
SRM's calculations and provide technical details on how 
to correctly enforce gauge-invariance for observables. Of 
course, thermal conductivity and any other linear re- 
sponse for a physical observable must be gauge-invariant 
independently of the means used to calculate it. In 
the framework of Kubo diagrammatic approach (both 
for a normal metal and a superconductor close to Tc), 
the relevant diagrams for thermoelectric response con- 
tain three sources of the vector potential/magnetic-field 
dependence: (A) Green functions/propagators, (B) heat 
current vertex, and (C) electric current vertex. The re- 
sulting expression is gauge invariant only if all these three 
sources are taken into account in a consistent fashion and 
within a specific gauge. In Ref [13] it is claimed that the 
contribution (B) cancels the contribution (A) calculated 
by Ussishkin in Ref [4]. However the consistency be- 
tween calculations and gauge choices in the two parts of 
the same physical quantity is not discussed. Most impor- 
tantly, contribution (C) is not even mentioned by Sergeev 
et al., which makes their conclusion erroneous. This con- 
trasts sharply with our approach, where we follow the 
work of Galitski and Larkin and from the outset 
use Landau representation that automatically takes into 
account all sorts of couplings to the magnetic field in a 
manifestly gauge-invariant fashion. 

We have found that the use of Landau basis greatly 
simplifies the calculations and automatically circumvents 
the need to keep track of gauge invariance at intermedi- 
ate steps. However, near T^, the use of this method is 
not required and correct result can be obtained via per- 
turbative expansions of all elements of the diagram in 
the vector potential. In this approach, care is required to 
keep the gauge invariance under control at all stages of 
the calculation. The correct starting point is to consider 
the correlator 



n„Mk,Q) = (j?(k + Q)jK-k)), 



(1) 



where Q is the momentum associated with the vector 
potential: A(r) — ae'*^ "". Setting k = implies that the 
momentum Q comes into the heat vertex, while setting 
k = — Q corresponds to Q coming out of the electric 
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vertex. The limit k — and Q needs to be taken 
in the end. 

Now we expand all terms (two vertices and two prop- 
agators) to linear order in A. In the most general 
form, this procedure generates four contributions II^™) 
(m = 1, . . . , 4) of the form (we retained the leading order 
in k and Q): 



tt(™) - c(™) I. a 



r 



(m) 



(2) 



Since 5^™^^ ^ 0, the relative contribution of these four 
terms retains the value of k. This is a reason why it is 
important to specify k when calculating each of the four 

gauge-noninvariant pieces. Of course, J2m ^a^jfj. ~ 
dicating that there is a well-defined limit nQ^(0, Q). The 
contributions T^^^^ are generally non-zero and explic- 
itly gauge-dependent, but their sum is certainly gauge- 
invariant: 



(3) 



where — L{ien + ii^iy,<z) and L_ = L{ien,q) are the 
fluctuation propagators. After summing over the Mat- 
subara frequencies, e„, and performing an analytic con- 
tinuation in uji, we come to the following result for the 
off-diagonal thermomagnetic response: 



/3, 



xy 
C.p. 



dL,LO^^\L{LO,q)flmLiLO,q), (4) 



where B{io) — coth(w/2T') is the bosonic equilibrium dis- 
tribution function, and 



L{uj,q) = 



ST,Jttv 



(8/7r)(T-Tc)-ia; + Dg2 



(5) 



is the retarded propagator for Cooper pairs. Evaluating 
the integrals in Eq. (jH), one immediately finds a finite 
value for Nernst coefficient oc (T — Tc)''/^"^ even in the 
absence of the particle-hole asymmet_ry, in accordance 
with all previous works in the field [2|- 



Thermomagnetic response: bosons vs. fermions 

To clarify the relation between the thermomagnetic re- 
sponses of fiuctuating Cooper pairs in the vicinity of 
[Eq. dU] and that of normal electrons, [i^^ , we find it 
instructive to present the latter in a similar form: 

/3^^oc I d^pp^J dee^^\G{e,p)\HmG{e,p), (6) 

where F(e) = tanh(£/2r) is the fermionic equilibrium 
distribution function, and 



G{e,p) 



1 



e-^p + i/2t 



(7) 



is the retarded electron Green function. The similarly 
looking expressions (|4]) and (O lead to essentially dif- 
ferent results for the Nernst coefficient for electrons 
and Cooper pairs. For free fermions in the absence of 
the particle-hole asymmetry, momentum integration in 
Eq. ^ gives an e- independent constant, and the subse- 
quent energy integral vanishes by oddness. On the other 
hand, for fluctuating Cooper pairs the integrals in Eq. ^ 
lead to a finite Nernst coefficient even in the absence of 
the particle-hole asymmetry. The same conclusion can 
be readily achieved within the time-dependent Ginzburg- 
Landau approach Q. 



Role of magnetization 

Sergeev et al. [H claim that the contribution of mag- 
netization is not relevant for calculating the Nernst co- 
efficient. This statement evidently contradicts the well 
established theory of thermomagnetic effects (see, e.g., 
[l5|). According to the theory, the experimentally mea- 
sured value of the thermoelectric tensor Pap is given by 
a sum of the kinetic contribution, /3q/3, (which can be 
calculated using the Kubo-like approach, see above) and 
a thermodynamic contribution due to magnetization M: 



(8) 



The importance of the magnetization heat current for 
the Nernst effect was first demonstrated by Obraztsov 
in 1965 16|. Later, the significance of this contribu- 



tion to thermomagnetic response of electron systems has 
been emphasized by a number of authors in relation to 
the integer quantized Hall effect, interacting electron gas 
in a quantizing magnetic field [l7|, and the fluctuation 
Cooper pairs contribution to the Nernst effect within 
Ginzburg-Landau formalism [^, The crucial role of 
the magnetization contribution to the Nernst effect has 
recently been demonstrated in Refs. 0, S] where it can- 
cels the otherwise divergent kinetic contribution /3xy at 
low temperatures above ifc2, thus being eventually re- 
sponsible for the implementation of the third law of ther- 
modynamics. 

SRM argued that magnetization effects do not con- 
tribute to the heat transport in a magnetic field, and the 
second term in Eq. ([5]) should be omitted. This statement 
contradicts all known theories and is incorrect. As a side 
remark, we mention here that even if the magnetization 
contribution is omitted, the Nernst effect in fluctuating 
superconductors would still be giant provided that the 
error in calculating made in Ref. [12[ is corrected. 

Finally, it also should be mentioned that the split- 
ting of the Nernst coefficient ([8]) in a kinetic (Kubo- 
like) and thermodynamic (magnetization) contribution 
is completely analogous to that in the Hall conductivity: 



'xy 



xy 



18|. 
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Summary 

To summarize this part: pretending to have developed 
a gauge-invariant microscopic approach to heat trans- 
fer, the article starts with a manifestly gauge non- 
invariant expression for the heat current [Eq. (5)] even 
for free electrons. It is precisely the incomplete account 
for the gauge invariance for why the authors of Ref. 
came to the wrong conclusion that the nonzero Nernst ef- 
fect due to fluctuation Cooper pairs necessarily requires 
presence of the particle- hole asymmetry in the underlying 
electronic band structure in a metallic phase. 



ON THE QUALITATIVE FORMULA FOR THE 
NERNST COEFFICIENT 



We take the opportunity here to discuss in more detail 
the qualitative expression for the Nernst coefficient, 



N = 



nq-^c 



dii 



(9) 



which was suggested in our Letter [3] and that was also 
criticized by SRM ^. 

We note that the main purpose of our work was to 
present a technically challenging microscopic calculation 
of the Nernst coefficient for an extended range of mag- 
netic fields and temperatures. These microscopic calcula- 
tions do not rely on any phenomenology such as Eq. ^ , 
but are strongly supported by it. The main purpose of 
including Eq. ([9]) in our paper was to provide a simple 
intuitive argument behind the giant Nernst effect in a 
fluctuating superconductor, which may be of value to the 
wide audience of PRL. And we believe that phenomeno- 
logical Eq. ([9]) is a new interesting (qualitative) result 
that does convey the desired message. 



Application to fluctuating Cooper pairs 

Arguing for the validity of phenomenological Eq. ([9]) 
we used the notion of a drift velocity and ignored im- 
purity scattering of the charge carriers, keeping in mind 
that the fluctuating Cooper pairs are not scattered by 
elastic impurities (all information about such scattering 
is included in the value of the effective coherence length 
or equivalently in the effective mass of the Cooper pairs). 
Equation Q clarifies the physical origin of the observed 
anomaly and, moreover, sheds light on where to look for 
the giant Nernst effect: in systems with chemical poten- 
tial strongly dependent on temperature. 

Now we clarify the issue related to the chemical poten- 
tial of fluctuating Cooper pairs, /ic.p.(r), which we intro- 
duced as an auxiliary concept within the phenomenolog- 
ical derivation of the Nernst effect. Indeed, it is known 



that in the thermodynamic equilibrium, the chemical po- 
tential of a system with a variable number of particles is 
zero, with photon and phonon gases being the textbook 
examples. A naive application of this "theorem" to fluc- 
tuating Cooper pairs "gas" has lead SRM to a wrong 
conclusion that //c.p. =0 [l]. However, a delicate issue 
concerning Cooper pairs is that they do not form an iso- 
lated system but are composed of fermion quasiparticles 
which constitute another subsystem under consideration. 
According to the same textbook discussion [i^j , in a mul- 
ticomponent system, the chemical potential of the ?'th 
component, /i^, can be deflned as the derivative of the 
thermodynamic potential with respect to the number of 
particles of z-th sort: 



py.Ni ' 



(10) 



provided the numbers of particles of all other species are 
fixed, Nj^i — const. In deriving the condition for ther- 
modynamic equilibrium, one should now take into ac- 
count that creation of a Cooper pair must be accompa- 
nied by removing two quasiparticles from the fermionic 
subsystem. This leads to /ic.p. ~ 2/i„ = 0, where /i„ is 
the chemical potential of quasiparticles. Therefore, the 
equilibrium condition does not restrict /ic.p. to zero, even 
though the number of Cooper pairs is not conserved. 

The value of /ic.p.(T) = T - |7| can be found from 
Eq. pH)) using the explicit temperature dependencies of 
fluctuation part of Gibbs potential flfi and concentration 
of fluctuating Cooper pairs. After this Eq. ([9|) reproduces 
[3| the known results for fluctuating Nernst response near 
the classical superconducting transition point. This style 
of reasoning is in close analogy with the "pedestrian" 
approach to fluctuating paraconductivity [6], which em- 
ploys the standard Drude formula with the assumption 
that the Ginzburg-Landau relaxation time plays the role 
of the scattering time. In this sense, Eq. (|9]) is an ana- 
logue of the classical Drude formula for electrical con- 
ductivity. This analogy also suggests that the domain of 
validity of Eq. ([9|) is constrained to temperatures close 
to Tc- In the opposite low-temperature limit, only micro- 
scopic calculations (such as presented by us in Ref. .7:]) 
are reliable. 

What concerns the applicability of Eq. ([S]) to a normal 
metal, one can easily check that it readily reproduces the 
Sondheimer's result. 



CONCLUSIONS 

To summarize, we list below certain points of the crit- 
icism of Ref. [13 together with our responses: 

1. SRM claim that "the linear response calculation 
of pxy does not require any magnetization correc- 
tion" and present their own interpretation of var- 
ious contribution to the heat current. — We cer- 
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tainly disagree with this statement. In such a del- 
icate and controversial issue as heat transport, we 
consider counterproductive to discuss interpreta- 
tions. Instead, one should either microscopically 
derive the expression for or check existing ones 
for possible inconsistencies. The crucial role of 
the magnetization contribution to Pxy bas recently 
been demonstrated in Refs. 0, S]: In particular, 
at T ^ it cancels the otherwise divergent Pxy, 
thus ensuring the implementation of the third law 
of thermodynamics. Therefore, omitting the mag- 
netization contribution to f5xy as suggested by SRM 
will inevitably violate the fundamental law of ther- 
modynamics. 

2. SRM claim that the particle- hole asymmetry of the 
single-electron spectrum is necessary for a nonzero 
Nernst effect due to fluctuating Cooper pairs. — 
This wr ong statement is solely based on the results 
of Ref. [12| which contains multiple errors, as dis- 
cussed in details above. 

3. SRM claim that our results generalized to repul- 
sive interaction in the Cooper channel would yield 
I3xy significantly exceeding that for non-interacting 
electrons. Further they claim that, "certainly, this 
effect is not known" . — We believe that the first of 
these statements is correct and the Cooper-channel 
correction to Pxy is indeed nonzero in the absence 
of the particle-hole asymmetry. Similar effects have 
been theoretically predicted for fiuctuation dia mag - 



tal importance, are completely unfounded, because the 
underlying theoretical work of SRM contains 12| multi- 
ple errors and inconsistencies. For this reason, there is 
no need to revise the existing theories of thermoelectric 
response in superconductors and in particular results of 
a microscopic analysis presented in our Letter Q; they 
remain intact. 



netism [2l| - |23l | and fluctuation conductivity [24 1. 
However we disagree with the second statement of 
SRM: The relative values of various contributions 
to Pxy should be considered individually for each 
experiment. 



4. SRM claim that, "according to textbooks [l9|, V/i 
should always be included in the effective electric 
field". — That is true, the textbook condition of 
vanishing current, _E — V/i/e — 0, is precisely our 
initial assumption of electroneutrality employed in 
the phenomenological treatment. 

5. SRM claim that "a thermodynamic value of /^c.p. is 
always zero, because a number of pairs is not con- 
served" . — This statement is wrong: The chemical 
potential of fluctuating Cooper pairs is nonzero in 
equilibrium since they do not form an isolated sys- 
tem. 
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